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ji ■ 1 Introduction 



We consider the Sturm-Liouville equation 

y" (x) + {X - q(x)}y(x) =0 (0 < x < oo) (1.1) 



^ ' with a boundary condition 

^t" ! y(0)cosa + y'(0)sina = 0, (1.2) 

^ ■ 

■ A being the complex spectral parameter. As usual, a is real and the potential q is real- valued and locally integrable 

' on [0, oo). We further assume throught the paper that q decays as x — ► oo in the sense that 

o 

O- g eL(0,oo). (1.3) 

+-> 

03 



Let us write A = z 2 , where < arg z < ir when < arg A < 2ir. Then ( [l. 3| ) implies that there is a solution 
tp(x, z) of ( ]1. l[ ) such that 

tp(x, z) ~ exp(zza;), ^ (s, z) ~ exp(izx) (1. 4) 

as x — > o o, and t/;(a;, z) is analytic in 2 for im z > J?], Theorem 1.9.1]. Then ^(x, z) is the Weyl L 2 (0 , oo) solution 
of ( |l. l|) w hen A is non-real and it forms the basis of the Weyl-Titchmarsh spectral theory of ( 1. 1 ) || Chapter 
9] 7 |17|,|18|. A central result of this spectral theory is the existence of a spectral function p a (jJ,) (— oo < \i < oo) 
which is piecewise constant in (— oo, 0) and locally absolutely continuous in [0, oo) with p a (n) > \ \L7[ section 5.7], 
|^8|, p. 264]. In particular, (1. 4) leads to the Kodaira formula 



L8|, p. 264]. In particular, ( 1. 4 ) leads to the Kodaira formula 

KPail 1 ) = M 1/2 / I *(A* 1/2 ) I 2 (M>0) (1.5) 

where 

^(z) = ip(Q, z) cosa + tp (0, z) sina (1-6) 

fl3[ p. 940]. Since the only possible eigenvalues of the problem ( ]l. l| )-( |l. 3| ) lie on the negative real A-axis, #(z) 
has no zeros for < arg z < tt expect possibly when arg z — ^tt. 
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the Weyl-Titchmarsh function m a (X) @, Chapter 9],jl7], Chapter 2] also involves 'J in the 



In addition to 
form 

m a {\) = {ip'(0, z)cosa--0(O, z) sina}/*(z) (1. 7) 

again with < arg z < n. Now m a (A) is related to the Green's function and to the resolvent operator of (|l] 
1 D-(B-- 2j) in the Hilbert space L 2 (0, oo), and the question arises whether these three spectral objects have analytic 
continuations into the so-called unphysical sheet tt < arg z < 2ir. As far as the Green's function and resolvent 



are concerned, this question can be posed, not only for (L. !)-( 



L. 2), but also for the corresponding Schrodinger 



equation in two or more dimensions. However, in the case of (1. 1) itself, it is a question of the analytic continuation 
of ip(0, z) and t/>'(0, z) in fll~6j) and §7% . 



Analytic continuation into the strip ~|a < im z < was established in || subject to a strengthening of ( |l. 3| ) 



to 



q(x) = 0{e- ax ) {x -> oo) 



(1. 8) 



for some a > (see also JlJ, section 2.2 
Allowing a to be arbitrarily large in ( 



1. 



and we refer again to || for a description of earlier work in this direction, 
leads to the class of super-exponentially decaying potentials for which 



with some f(x) — ► 
specalisations of (1 



q{x) = 0{e-"t^) (x -> oo) 

and then we have analytic continuation into the whole of im z < 0. JTof , 
where again there is analytic continuation into the whole of im z < are 



(1. 9) 



121. Two other 



q(x) 



with p(x) periodic and 



q(x) — (const. )x N e 



where jV(> 1) is an integer Q. 

Once analytic continuation has been effected, the possibility is opened up of *&(z) having zeros in the unphysical 
sheet. Such zeros are called resonances and, by < \l. 7\) , they a re sin gular spectral points associated with the Green's 
function and resolvent operator. For potentials of the class (|l. 9|) , the asymptotic distribution of resonances was 
obtained i n [flo[ and, by another method, also in (l2) along with other results on the location of resonances. In 
particular |12, Theorem 3.8], there is a resonance-free strip —b < im z < subject to q having a suitably small 
norm. 

All these existing results require exponential decay of the potential q. In this paper, we allow q to have only 
power decay 0{x~ 1 ) (x — > oo) for some 7 > 1 and, under further conditions on the analyticity of q, we establish 
analytic continuation of *&(z) into a sector 2tt — 6q < arg z < 2tt of the unphysical sheet. The necessary construction 
is given in section 2. Then in sections 3 and 4 we show that our methods lead to certain resonance- free regions which 
are adjacent to part of the real z-axis. Finally, in section 5, we discuss the numerical computation of resonances 
lying in the complement of our resonance- free regions. 



2 Analytic continuation 

The method which we develop in this section for continuing tp(x, z) and tp (x, z) analytically into im z < is based 
on the integral equation by means of which ( 1. 4 ) is proved [0, sections 1.3 and 1.9]. Thus we begin by writing (|l] 
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w = ( ! 7/: ) ( i i ■ p. i) 



1) (with A = z 2 ^ and y = ip) as a first-order system in a standard way by defining 

i -te, ( i -*/* V ^' 

2 I 1 i/* J I V> 
Then 

<*-'={( ° -L) +e ("i 1 "i 1 )}'*'' < 2 ' 2 > 

where 

Q = ~i<7/z, (2. 3) 

and the corresponding integral equation is 

W(x,z)=e 1 + / Q(t)K(t - x, z)W(t, z)dt, (2.4) 



where 

ei = f J, #(s,*)=( _ £ 2isz _ e 2isz )• < 2 - 



Iteration of (2. 4) gives 

oo 

W(a?,«) =e 1 +J2W n {x,z), (2.6) 
i 

where 

/>00 

W n {x,z)= Q{t)K(t-x,z)W n -i(t,z)dt (2.7) 



and Wo (a;, z) = e±, provided of course that the infinite integral converges. We note that, in terms of the components 
of W n and W n - U (O is 

( U v n n \ (x,z) = J™ Q(f)K-i (*>*)+«»-! («,«)} ( _ &2 ^ x)z ) dt. (2. 8) 

Also, the transformation ( 2. 1 ) back to ( |l. l|) gives 

^ = e lxz (l + f> n (x, z) + v n (x, z)}^ , (2. 9) 

i>' = ize lxz (l + f2{u n (x, z) - v n {x, z)}^j . (2. 10) 

In what follows, we write 

\W n \ = \u n \ + \v n \. (2. 11) 

We now introduce the more detailed conditions on q that we require. 
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Condition 2.1 We suppose that the real-valued function q(x) can be extended into a sector S of the complex 
£— plane as an analytic function as follows. 

1. is regular in a sector S defined by —9q < arg £ < 9q and £ ^ 0, with some 9q such that < 6q < n. 

2. There are constants 7 (> 1) and k such that 

I q(0 \<k\s r (2. 12) 

as I £ |— > 00 and £ 6 5. 

Our met hod i nvolves extending also the definition of W(x, z) into the complex £-plane. To do this, we write 
t = x + s in (|2. 7j ) and consider the iterative definition 



Q(£ + s)K(s, z)W n -i(£ + s,z)ds 



(2. 13) 



for £ € 5, with Wo(£ ; z) = ei. In the following l emma we give a simple estimate for the size of W n in order to deal 
with the convergence of the infinite integral in ( |2. 13| ). 



Lemma 2.2 Let q satisfy Condition 2.1. Let (eS and let im z > in l[2. dj ) and \2. IS ). Then for n>0 

1 



(2. 14) 



and W / Tl (£, z) is a regular function of ^ in S. 



Proof. We note that the infinite integral in ( 2. 14 ) converges beca use of ( 2. 12 ). The lem ma is clearly true w hen 
n = and, proceeding by induction on n, we use the form of ( 2. 8 ) which corresponds to ( 2. 13 ). By ( |2. ll[ ) and 
(|2. 14|) (with n - 1), this gives 



1^(^)1 < ^riQ^^lj^TyXlo 



(n 



s + cr) I <i<7 

Tl-l 



from which (2. 14) follows. 

To deal with the regularity of the W n , we note that ( [2. 12[ ) and ( [2. 14[ ) imply that the infinite integral in (^] 
[l3| ) converges uniformly with respect to £ in any closed bounded region Si C S. Thus the regularity in S of W n 
follows from that of W n -i, and the lemma is proved. 

The next step is to re- write ( [2. 13 ) in a form which does not require im z > and which therefore provides the 
analytic continuation of W n (£, z) (as a function of z) into the lower half of the z— plane. At this stage we restrict 
£ so that re £ > 0, the reason being given in the proof of the following theorem. Ultimately we specialise £ to be 
the positive real variable x. 
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Theorem 2.3 Let q satisfy Condition 2.1. Then, for all £ and z in S with re £ > 0, 

i r°° s s 

W n (£,z) = - Q(£+-)K(8,l)W n -i(Z+-,z)ds, (2.15) 

Z Jq z z 

and the series 

oo 

W%z) =ex + J2 W n& z ) ( 2 - 16 ) 

1 

defines a regular function of z in S which, when £ = x, continues to satisfy the differential equation 1 ). 
Proof. We suppose first that im z > 0, so that ( b! 13 ) holds. We consider the contour integral 

Q(£ + iJw^-iK + -> *)*n (2- 17) 

z z 



c 



where C is the closed contour in the complex plane formed by the positive real axis, the line through z from to 
oo, and the smaller part of the circle | r\ |= R. The assumption that re £ > guarantees that the point £ + rj/z 
lies in S, and therefore the integrand in ( |2. 17| ) is defined as a regular function of r\ within and on C. Then, by 
Cauchy's Theorem, the value of ( |2. 17| ) is zero. Thus (2. 15) follows from ( [2. 13 ) when R — > oo, provided that the 
contribution to (2. 17) from In |= R tends to zero. 

By (2. 3), (2. 5) and ( |2. 14| ), this contribution does not exceed in modulus 



(const. )R 



v 

a; 
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2; 



n-l 



(2. 18) 



in which r\ = i?e 4e . Since < arg z < 6*o (< 7r), it is easy to check that 

\X + iY + ^\ 2 > I cosf? |)(^ 2 + i? 2 / U I 2 ) - ^ 

z 2 

where X = re £ + s, Y = im £ and k = (1 + 3 cos 2 flp )/■ [(!— | cos#o |)(1 + 3 | cos0 o |)}. Then, since 7 > 1, the 
9— integrand in ( 2. 18| ) is 0{R i' _ (" _1 )(t _1 )) ; and hence (2. 18) tends to zero as R — * 00 for all n > 1. This proves 
Q2. 15|) for im z > 0. 

We turn now to im z < and we show that ( |2. 15| ) continues to provide an iterative definition of the W n (£, z) 
as regular functions of z. We note that, when re £ > and im z < 0, the point £ + s/z in ( [2. 15j ) continues to lie 
in S. An induction argument similar to that used for ( |2. 14 ) shows that 



1/0 poo \ 71 

W n (Z,z)\<-l^—j o I Q(£ + -)|<fe) . 



(2. 19) 



Again, as for ( [2. 13|), the infinite integral in ( 2. 15 ) converges uniformly with respect t o z in any closed bounded 
region Si C S, by ( [2. 12 ). Hence e ach W n {j, z) is a regular function of z in 5. Further, ( 2. 19 ) also guarantees the 
uniform convergence of the series ( 2. 16 ) with respect to z in S\, and hence W(£, z) is also a regular function of z 
in S. 
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Finally, we show that W(£, z) satisfies ( 2. 2| ) in the more general form with £ in place of x. In (2. 15), we sum 
for n going from 1 to oo and we write s = z(t — £) to obtain 



W(t z) = ex + J Q(t)K(z(t - 0, l)W(t, z)dt, 



(2. 20) 



where oo denotes the point at infinity on the line through £ in the direction of the vector 1/z. The interchange 
of integration and summation involved in ( 2. 20| ) is justified by means of 02. 19D . Differentiation of (|2. 2C|) with 
respect to £ now recovers (2. 2) with £ in place of a;, and the proof of the theorem is complete. 
We note that, in terms of the components of W n and W n -\, (2. 15) is 



(2. 21) 



corresponding to ( |2. 8| ), and we have used ( J2. 3| ). Here (2. 21) is valid for all £ and z in 5 with re £ > 0, and then 
( |2. 9| ) and ( |2. 10| ) provide the desired analytic continuation of ip(x, z) and ip (x, z) into the lower half of the sector 
S. 



3 Resonance-free regions 



The basic result on non-resonance which follows from ( 2. 9 ), ( 2. 10 ) and ( 2. 21 ) is given in the next theorem. 
Theorem 3.1 Let z £ S with im z < and z ^ icota. Let 



|g( S /z)|d S <|z| 2 log(l + ( 5- 1 ), 



(0 < a < tt/2 



where 

5= { 1 

1 | cos a — iz sin a | / | cos a + iz sin a | (tt/2 < a < it) 

Then ^(0, z) ^ and z is not a resonance. 
Proof. By ( p~6| ), ( p7~9| ) and ( p. 10| ), we have 



(3. 1) 
(3. 2) 



^(z) = (cosa + izsina) I 1 + ^^{u„(0, 2) + Zv n (0, z)} 



(3. 3) 



where Z = (cosa — izsina)/(cosa! + izsina). It i s eas y to check that, since im z < 0, | Z |< 1 for < a < tt/2 
and I Z |> 1 for tt/2 < a < ir. Hence, with S as in ( 3. 2 ), 



£{«n(0,*) + Zl> n (0,*)}| < ^(|u„(0, z) | + | w n (0, z) |) 



; M <'*!>( I - I - * / I g(s/z) I da J - 1 
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by (2. 19) and (2. 21). It now follows from (3. 3) that $f(z) is non-zero if z ^ icota and 



cxp (j z | 2 jf 



I q( s / z ) \ ds) —1 < <5 _1 , 



and the latter is guaranteed by ( |3. l[) . 



Let us note that, with the change of variable s =\ z \ t, ( 3. 1 ) can be written as 

\q(te lS ) |^<|z|log(l + ,r 1 ) 



(3. 4) 



where — 2tt — arg z > 0. The condition (3. 4) defines a region of the complex plane within which there are no 
resonances, and the nature of this resonance-free region depends on the nature of q. Befo re we turn to detailed 
examples, we give one general property of resonance-free regions which is a consequence of ( 3. 4 ). 

Corollary 3.2 There are real numbers R\ (> 0) and 9\ (0 < &\ < n) such that the sectorial region \ z \> Ri, 
2-7T — B\ < arg z < 2tt is resonance-free. 



(3. 5) 



(3. 6) 



Proof. Suppose first that < a < 7r/2, so that 5 = 1 in ( 3. 2 ). We choose Ri so that 

poo 

fll>(log2)- 1 / \q(t)\dt. 
Jo 

Then, by continuity in 8, we have 

/>oo 

/ I q(te i9 ) I dt < i?ilog2 

for 9 in some range (O,0i) with 0i > 0. Hence ( [Til) holds for | z \ > R±, and the corollary is proved for this range 
of a. 

Next suppose that ir/2 < a < it. Then, with S as in ( |3. 2| ) and cot a < 0, it is easy to check that 

<5 < (1 + sin0)/ | cos0 | 

where again arg z — 2ir — 6. We choose R\ as in (3. 5) but, in place of ( 3. 6| ), we can say that 

cos 



I" 3 | q(te ie ) \dt<R x log ( 1 + J mS0 [ ) < R, log(l + 5- 
Jo V l + sm0y 



for in some range (O,0i) with 6\ > 0. Hence (3. 4) again holds for | z |> as required. 

Corollary 3.2 provides theoretical support for an observation by Aslanyan and Davies || concerning the numer- 
ical computation of resonances for the potential 



q(x) = x 2 exp (— ex 2 ) , 



(3. 7) 



where e (> 0) is a small parameter. In p. 16 and Table 10] it is noted that there are resonances very close to 
the positive real axis but, at a certain point, they turn sharply away into the lower half plane. Now (3. 7) satisfies 
Condition 2.1 with 0q < 7r/4 (see also Example 4.6 below), and the existence of the sectorial region in Corollary 3.2 
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precludes as a general feature the occurrence of resonances close to the positive real axis beyond a certain distance 
from the origin. 

Corollary 3.2 can also be related to J8|, Theorem 1] concerning the localization of spectral concentration points to 
a bounded interval on the real spectral axis (see also || section 2]). Insofar as spectral concentration is associated 
with resonances located near to the real axis, Corollary 3.2 provides anot her p roof that spectral concentration 
points are confined to a bounded interval for a class of potentials satisfying (1. 3). 

In the Dirichlet case a = of ( |l. 2| ), there are additional non-resonance results like Theorem 3.1 and Corollary 
3.2 but with, in the corollary, the vertex of the sector at the origin. 

Theorem 3.3 Let a = in \l. ^j ) and let q satisfy Condition 2.1 with 

7 > 2 (3. 8) 



in (2. li ). Let z 6 S with im z < 0, and let 

/•OO 

/ s | q{s/z) \ds<\z | 2 log 2. 
Jo 

Then ip(0, z) ^ and z is not a resonance. 



(3. 9) 



Proof . We note that ( |3. q ) guarantees the convergence of the integral in ( 3. 9 ). In ( 2. 21 ), we use the inequality 



I < 2s to obtain 



(u n + v n )(£, z) \<\ z [ 



s I <?(£+-) II Ki-i +«ti-i)(£+ -, z ) I ds. 



Then, as for (2. IE), an induction argument gives 



I (U n +V n )(€,z) \< — 
n 



This inequality is used in ( p. 9]) and (1. 6) (with a = 0), and the theorem follows from (3. 9) in the same way as 

Theorem 3 .1 fo llowed from ( |3. l| ). 

As for ( 3. 4 ), the change of variable s =| z \ t in (|3. 9| ) leads to 



(3. 10) 



/ t I q{te w ) | dt < log 2 
Jo 

and this in turn leads immediately to the next corollary. 
Corollary 3.4 Let a = in ( fi. ^ ) and, in addition to (3. t ), let 

/>oo 

/ t I q(t) | dt < log 2. (3. 11) 

Jo 

Then there is a real number 9\ (0 < 6\ < ir) such that the sector \ z |> 0, 2ir — B\ < arg z < 2n is resonance-free. 
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The condition fl3. llj) can be related to the condition 



/ t | q(t) | dt < 0.1735 (3. 12) 

Jo 

H (2.16)] which is shown in Q (by quite different methods) to imply the absence of any spectral concentration 
points on the positive spectral axis (0, oo). The smaller the value of the integral in ( |3. 11| ), the larger 9\ can be, 
and the further away from the real axis are any resonances pushed. Thus, in the case of ( |3. 12| ), any resonances 
are too far from the real axis to produce spectral concentration fi section 3(iv)]. 



4 Examples 

We consider now some examples of q which sho w in m ore detail the type of region that arises from ( |3. 4| ). We 
keep to the case < a < vr/2 for which 5 = 1 in fTJ): in the other case, S — > 1 as | z |— > oo and the regions are 
asymptotically similar for large | z \. 

4.1 Example q(x) = c(x + a)~ 7 

where 7 > 1, c and a are real and a > 0. In Condition 2.1, we take g(£) = c(£ + a) -7 with, if 7 is not an integer, a 



cut in the £— plane from —a to — 00 along the real axis. Thus we can take 0q = n. The integral in (3. 4) is now 







\t + ae- w dt=\c\ i (t 2 + 2at cos 9 + a 2 )-~> /2 dt = 1(6) (4.1) 



say. Hence 1(9) increases from | c | a 7+1 /(7 — 1) to 00 as 9 increases from to it, and (3. 4) becomes 



\z |>J(0)/log2. (4. 2) 

Thus we have a resonance-free region which lies in the lower half of the the complex plane, bounded by a curve 
which starts at the point | c | a 1+1 / (7 — 1) on the real axis and recedes from the origin as 9 (= — arg z) increases 
from to 7T. The region is of course on the side of the curve remote from the origin. When 7 = 2 in particular, the 
integration in ( |4. l| ) can be performed and (4. 2) becomes 



I ^ I > — j — (6» = -argz). 
a log 2 sin 9 

Thus the boundary curve in this case is asymptotic from above to the line im z — — \ c \ 7r/(alog2) as 9 — > n. 
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4.2 Example q(x) = c(x n + a 

where n (> 2) is an integer, 777 > 1, c and a are real and a > 0. This is similar to Example 3.1 but now 9q = ir/n. 
The integrand in ( [l. l|) is replaced by (t 2n + 2a n t n cosn9 + a 2 ™) -7 / 2 , and 1(9) increases to 00 as 9 — > ir/n. 
In the case when n = 2 and 7 = 2, we find that 

1(9) = (tt I c I /4a 3 ) sec 9, 
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and the resonance- free region ( 4. 2| ) is the quadrant 

re z > 7r | c | /(4a 3 log 2), im z < 0. 
Also in this case, the Dirichlet condition ( |3. 10 ) gives 

26>/sin2(9 < 2(a 2 / | c |)log2. 

Thus, on the assumption that | c |< 2a 2 log 2, this being (|3. ll[ ), the value of 6\ in Corollary 3.4 is the solution of 

20i/sin20 1 = 2(a 2 / | c |) log 2. 

4.3 Example q(x) = c{(x — w)(x — w)}~ 1 

where 27 > 1, w 7^ and < arg w < tt. This again is similar. Here 8q = arg w (= say ), and the integrand in 
( [4. 1| ) is replaced by 

{(t 2 - I u> | 2 ) 2 - 4 I w I t(t- I w |) 2 cos0cos6> + 4 I w | 2 i 2 (cos6» - cos^) 2 }" 7 / 2 . 

Again 1(8) — > 00 as 6 — + because we approach a singularity at t =| «; | in {...}~ 7 / 2 . However, 7(6*) is not 
necessarily monotonic unless cos0 < 0. 

We conclude this group of examples by noting that similar remarks apply when q is a product of terms already 
considered with differing values of a, w, and 7 and, indeed, when q is a ratio of two such products. We give one 
example of this more general type for future reference in Section 5. 

4.4 Example q(x) = c(x — l)/(x + l) 4 . 
Here 1(9) in ( |4. l| ) and (|~|) is replaced by 



1(8) =| c 



p OO 

/ |te ie -l I /(£ 2 + 2£cos(9 + l) 2 di. (4.3) 
Jo 



Now the boundary curve of the resonance- free region(4. 2) starts at the point 0.36 | c | on the real axis and, since 
1(9) ~ (const. )(sin6*)~ 3 when 8 — > tt, the curve is asymptotically like | re z | 2 = (const.) | im z | 3 (see also Figure 1 
below). 

Next, we turn to examples with exponential decay which are also covered by Condition 2.1. 

4.5 Example q(x) = 2e~ ax smx (a > 0). 
In Condition 2.1, we take 

g(£) = i ( e -( a +*)« - e -( a -^ . (4. 4) 

Since | e~^ a±l ^ \= e - Qre « ±im «, ( |2. 12| ) is certainly satisfied if 

8 a < tan -1 a. (4. 5) 
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By (4. 4), the left-hand side of (3. 4) does not exceed 



/>oo 

2 / exp {(— a cosO + sin#) t} dt = 2 (a cos# — sin 6>) _1 . 
Jo 



Hence (3. 4) holds if | z \ (acosO - sin#) > 2/ log 2, or 



a(re z) + (im z) > 2/ log 2. 



Since 9q can be arbitrarily near to tan -1 a in (4. 5), we therefore have a resonance- free region in the lower half of 
the complex plane lying to the right of the line through the point 2/(alog2) on the real axis and with gradient —a. 

We observe that independent support for this gradient —a is provided by the quite different analytic continuation 
method developed in j|, Prop. 2.1]. This latter method constructs the analytic continuation of ^(z) into the whole 
of im z < except for poles at the points 

— —(v + mai) (| v |< m, m = 1, 2, ...) , 

v being an integer. Thus the line through the origin with the same gradient —a delineates a pole-free region for 
^(z) within which $f(z) is regular. The methods in B do not however lead readily to resonance-free regions. 



4.6 Example q(x) = cx m exp(— x n ) 

where m and n are positive integers. In Condition 2.1, we take 



and (|2. 12j ) is certainly satisfied if 



«(0=crexp(-D, 
O < 7r/ 2n - 



(4. 6) 



Now the left-hand side of (3. 4) is 



c I / t m exp (-t n cos nO)dt=\c \ (cos nO) 
Jo 



— (rn+l) / n 



where I = J u m exp(—u n )du. Hence ( 3. 4 ) holds if 



z | (cosn6l)( m+1 )/" >| c | 7/log2. 



(4. 7) 



In ( [4. 6| ), 6*o can be arbitrarily near to n/2n and hence, in ( 4. 7 ), we can let 9 increase from to 7r/2n. Thus ( [4. 7| ) 
defines a region in the lower half plane whose boundary starts at the point | c | I / log 2 on the real axis and recedes 
to infinity as 9 — > ir/2n. 

A typical example of (4. 7) is when m — and n = 2, in which case the boundary is the part of the rectangular 



hyperbola X — Y = (| c | //log 2) which lies in the fourth quadrant of the (X,Y)— plane and z = X + iY. 
Again the resonance-free region lies on the side of the hyperbola remote from the origin. Independent support for 
the nature of this boundary is provided by the findings of Siedentop |la| and Froese [0. In JIB] (where c = — 1), 
the first few resonances found computationally are already near to, but below, the line arg z = —it /A (see also |^, 
Example 6.4]) while, in p3|, the resonances are shown to be asymptotically near to this same line. 
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5 Computational resonance-finding 



We turn now to the numerical computation of resonances for explicit q, such as those in section 4, which satisfy 
Condition 2.1. One possible direct method is to compute the u n and v n (1 < n < N) recursively in ( |2. 21| ) and 
substitute the results into ( 2. 9 ) and ( 2. 10 ), the infinite series being truncated a t N w ith an error term. Then a 
zero-finding algorithm would be applied to the resulting approximation to ^ in ( |l . 6| ) . A similar procedure was 
applied successfully to the formulae for u n and v n in |4 when q has exponential decay. However, in our present 



situation of power decay, it has proved difficult to use ( 2. 21 ) when N > 2, repeated integration being involved, 
and in addition the error term for N = 1 is not small. 

Instead, we have computed resonances by the method of complex scaling. We refer to Simon Jl6| for a discussion 
of this method in relation to resonances and to Agmon jjj for a recent definitive account in a very general setting. 
The method of complex scaling is closely associated with ( [2. 20 ) and (|T~9) and, in fact, our approach in section 2 
provides an independent justification of the validity of this method for ( p.. 1 ), as we now describe. 

The transformation of ( [2. 20 ) back to ip(£, z) via ( [2. l| ) (with £ in place of x) gives 



d^/df + {z 2 - q(t)}^ = 
corresponding to ( |l. l| ). With £ in polar form £ = rexp(i^) (0 < cf> < 8q), we therefore have 

e- 2l4, d 2 4>/dr 2 + {z 2 - qire^)}^ = 



(5. 1) 



and, by (|TJ) and < \2. 19[ ), 

| V(re # , z) |= [exp{-r | z | sm(</> + arg z)}]{1 + o(l)}, 
where o(l) refers to r — > oo. It follows that tp(re lc ^, z) is an L 2 (0,oo) solution of (fTT]) if 

2tt - <j> < arg z < 2tt. (5. 2) 

Thus the zeros of ^(z) in (|l. 6| ) provide the eigenvalues z 2 of ( ^3. l|) on < r < oo with the boundary condition 

■0(0, z) cos a + ip (0, z) sin a = 



at r = 0. Here ( |5. 1| ) is said to be obtained from ( |1. l| ) by complex scaling, the scaling factor being e 1 ^ ^ section 
5] |l6|, section 3]. 

We have therefore applied a computational eigenvalue finder to ( |5. l[ ) with a suitable value of cj>. This 
locates eigenvalues z 2 and hence resonances in the sector ( |5. 2| ). We focus the discussion of our computational 
findings now on Examples 4.1-4.4 since it is potentials with only power decay which are the main object of this 
paper. 

We consider first 

q(x) =c{x 2 + l)- 2 , (5.3) 

being the case n — 7 = 2, a = 1 of Example 4.2. Here 9q = ir/2 but, if cf> in ( |5. l| ) is close to tt/2, the code in 
^l| reports unreliable results due to the sharp (but non-singular) maximum of the scaled | q{re 1 ^) \ near to r = 1. 
Accordingly we have chosen </> = 1.5. We have found no resonances satisfying (|5. 2| ) within the disk | z |< 10 
when c has the range of values —1, —5, —10, —15, —20. This is certainly consistent with the resonance-free quadrant 
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re z > it | c | /(41og2), im z < in Example 4.2, but there remains the open question whether resonances occur 
elsewhere in im z < 0. 

A similar example, but with a higher singularity located nearer to the real axis in the complex plane, is 

q(x) = c(x 6 + I)" 20 . (5. 4) 

Despite this extra feature, this example also produces no spectral concentration and no resonances. Here 9q = n/6 
and we have c hosen <fi — 0. 5. Th e values of c investigated were — 1, —5, —15, —20, —25, —30. The reason for choosing 
c negative in is to give q(x) a negative minimum, a property which in exponentially decaying 



examples is often associated with spectral concentration and resonances [|l| Section 6]. 
Next we consider Example 4.4 

q(x) = c(x-l)/(x + l)\ 

this time with c > to give the negative minimum (at x — 0). Here of course 9q = n and we have chosen <fi = 3.0. 
There is one real point of spectral concentration when c = 35 located at A = 0.26, for which y/X = 0.51, and we 
have tracked the corresponding resonance for a range of values down to c = 0.5. The resonance broadly recedes 
from the real z— axis as c decreases, and we give a selection of these findings in Table 1. For small values of c, re z 



c 


rez 


im z 


35 


0.50 


-0.06 


25 


0.65 


-0.23 


15 


0.57 


-0.44 


10 


0.42 


-0.55 


4 


0.06 


-0.68 


3 


-0.05 


-0.68 


2 


-0.21 


-0.66 


1 


-0.45 


-0.58 


0.5 


-0.67 


-0.45 



Table 1: Resonance z 



appears to increase rapidly in the negative direction, but arg z becomes too close to tt for the code jllj] to produce 
reliable values. 

In order to gain an idea of how Table 1 relates to the resonance-free region given by ( p~2l ) and we note 



that 1(9) contains a factor | c |. Accordingly, we have applied a scaling factor | c | to both 1(9) and the values in 
Table 1. The result is Figure 1, in which the diamonds denote the scaled resonances from Table 1, and the dotted 
curve denotes the boundary curve scaled to c = 1. Figure 1 confirms the general nature of our theoretical result 



(Li) 



We also mention that there are two additional similar strings of resonances: when c = 10 for example, there 
are resonances at z = —1.27— 1.39i and z = —5.05 — 5.31i in addition to the value in Table 1. These additional 
resonances, however, lie further from the resonance- free region than the resonance-string shown in Figure 1. 

Finally, we have also considered the example 

q(x) =c(x-l)/(x i + l) 
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Figure 1: 



for which 0q = tt/4 and we have taken <f> — 0.75. There is one real point of spectral concentration when c = 11 
located at A = 0.15, for which y/\ = 0.39. We have tracked the corresponding resonance from z = 0.39 — 0.03t 
when c = 11 as far as z = 0.73 — 0.59i when c = 3.2. For smaller c, the code H] again flags unreliability, but there 



is a corresponding picture to Figure 1 to similarly confirm the theoretical result (3. 4) 
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